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Abstract 

Semi-empirical interatomic potentials have been developed for Al, a-Ti , and 7-TiAl within the 
embedded atom method (EAM) formalism by fitting to a large database of experimental as well as 
ab-initio data. The ab initio calculations were performed by the linearized augmented plane wave 
(LAPW) method within the density functional theory to obtain the equations of state for a number 
of crystal structures of the Ti-Al system. Some of the calculated LAPW energies were used for 
fitting the potentials while others for examining their quality. The potentials correctly predict the 
equilibrium crystal structures of the phases and accurately reproduce their basic lattice properties. 
The potentials are applied to calculate the energies of point defects, surfaces, and planar faults 
in the equilibrium structures. Unlike earlier EAM potentials for the Ti-Al system, the proposed 
potentials provide a reasonable description of the lattice thermal expansion, demonstrating their 
usefulness for molecular dynamics and Monte Carlo simulations at high temperatures. The en- 
ergy along the tetragonal deformation path (Bain transformation) in 7-TiAl calculated with the 
EAM potential is in a fairly good agreement with LAPW calculations. Equilibrium point defect 
concentrations in 7-TiAl are studied using the EAM potential. It is found that antisite defects 
strongly dominate over vacancies at all compositions around stoichiometry, indicating that 7-TiAl 
is an antisite disorder compound in agreement with experimental data. 

PACS numbers: 61. 72,62. 20.Dc,64.30.+t,65.40 
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I. INTRODUCTION 



In the recent years, intermetallics alloys based on the gamma titanium aluminide TiAl 
have been the subject of intense research due to their potential for applications in the 
aerospace and automobile industries.— Such alloys have an excellent oxidation and cor- 
rosion resistance which, combined with good strength retention ability and low density, 
make them very advanced high temperature materials. A study of fundamental properties 
such as the nature of interatomic bonding, stability of crystal structures, elastic proper- 
ties, dislocations, grain boundaries, interfaces, as well as point defects and diffusion are 
therefore warranted in order to gain more insight into the behavior of these intermetallic 
alloys under high temperatures and mechanical loads. Over the past few years, numer- 
ous investigations, both experimental and theoretical, have been devoted to the study of 
such properties . 5l6l7l8l9ll ° l11112113 Accurate ab-initio studies of the structural stability, elastic 
properties, and the nature of interatomic bonding have been reported for 7-TiAl as well 
as other stoichiometric alloys of the Ti-Al system . 14 i 15 i 16 However, the application of ab- 
initio methods to atomistic studies of diffusion, deformation, and fracture are limited due 
to the prohibitively large computational resources required for modeling point defects, dis- 
locations, grain boundaries, and fracture cracks. Such simulations require large simulation 
cells and computationally demanding techniques such as molecular dynamics and Monte 
Carlo. Semiempirical methods employing model potentials constructed by the embedded 
atom method (EAM) 17 ' 18 ' 19 or the equivalent Finnis-Sinclair (FS) method 2 ^ are particularly 
suitable for this purpose. These methods provide a way of modeling atomic interactions 
in metallic systems in an approximate manner allowing fast simulations of large systems. 
Several studies applying these methods to a variety of properties of 7-TiAl such as pla- 
nar faults, dislocations, and point defects have been reported in the literature . 12121122 The 
effectiveness of semiempirical methods obviously depends upon the quality of the model po- 
tentials employed. Recent studie o 2 ^ 2412 ^ 2 ^^! 2 ?! 2 ?!??!? 1 !? 2 have shown that the incorporation 
of ab-initio data during the fitting of interatomic potentials can significantly enhance their 
ability to mimic interatomic interactions. For example, Bakes et alM examined the range of 
interatomic forces in aluminum using model potentials and ab-initio methods. They found 
that potentials that included ab-initio data during the fitting procedure could reproduce 
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ab-initio forces much more accurately than potentials fit to experimental data only. 

In the present work we explore the possibility of constructing a reliable interatomic po- 
tential for the Ti-Al system. To this end, we develop EAM-type interatomic potentials for 
7-TiAl and the component elements Ti and Al by fitting to a large database of experimen- 
tal properties and ab-initio structural energies of these phases. The ab-initio database has 
been generated by density functional calculations using the linearized augmented plane wave 
(LAPW) method within the generalized gradient approximation (GGA) for the exchange- 
correlation effects. The ab-initio data are used in the form of energy- volume relations (equa- 
tions of state) of various structures of Al, Ti, and 7-TiAl . The energy along the Bain 
transformation path between the Ll and B2 structures of TiAl has also been calculated in 
this work. 

While many EAM-type potentials have been reported for Al , 23 ' 28 ' 32 ' 33 ' 34 relatively few 
attempts have been made to create such potentials for T i 35 i 36 i 37 and TiAl . 21 i 22 i 38 i 39 i 40 i 41 i 42 
Titanium, like other transition metals, cannot be expected to follow the EAM model as 
accurately as noble metals usually do. In contrast, Al is known to lend itself to the EAM 
description quite readil y 23 i 28 i 32 i 33 i 34 The 7-TiAl compound is probably a borderline case. 
Some of the previous EAM potentials for 7-TiAl had a reasonable success in modeling 
lattice properties and extended lattice defects . 21 ^ 22 ! 38 On the other hand, Paidar et a/.— 
calculated deformation paths between different structures of TiAl by an ab-initio method 
and with a Finnis-Sinclair potential and found only a qualitative agreement between the 
two calculation methods. In this work we are trying to gain a better understanding of 
the applicability range of the EAM model for 7-TiAl . In particular, we want explore the 
possibility of obtaining a high-quality EAM potential for 7-TiAl by incorporating ab-initio 
data into the fitting database. Our results suggest that the quality of an EAM potential 
can indeed be improved this way, and that the potential proposed here can be useful for 
atomistic simulations of 7-TiAl . 
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II. CONSTRUCTION OF EAM POTENTIALS 



A. Database of fitted physical properties 

We have first constructed EAM potentials for pure Al and Ti followed by fitting a cross- 
interaction potential Ti-Al. The database of physical properties employed in the fitting 
procedure consisted of two categories. The first category is comprised of experimental data 
for the lattice constant, c/a ratio, cohesive energy, and elastic constants. For pure Al 
and Ti it also included the vacancy formation energy and the linear thermal expansion 
factors at several temperatures. The second set of properties consisted of ab-initio energy 
differences between various crystal structures. Such differences are necessary to ensure the 
correct stability of the experimentally observed ground state structures against other possible 
structures and to sample a large area of configuration space away from equilibrium. 

The ab-initio database consists of energy versus volume (EV) relations for various crystal 
structures. For Al, EV curves were computed here for the face centered cubic (fee), hexag- 
onal closed packed (hep), body centered cubic (bec), simple cubic (sc), diamond, and the 
110 structures. The Ll structure of Al is a defected fee lattice with a vacancy in the corner 
of each cubic unit cell. In the case of titanium, EV curves were generated in Ref. |44j for the 
hep, fee, bec, sc, and the omega (C32) structures. For the intermetallic compound TiAl, EV 
relations were obtained in this work for three structures: 110 (CuAu prototype), B2 (CsCl 
prototype), and Bl (NaCl prototype). Each EV curve typically consists of total energies 
for about 20-30 different volumes around the equilibrium volume. The c/a ratio of the 110 
structure has been optimized at each volume. As the ab-initio and EAM energies are at dif- 
ferent scales, all ab-initio energies for a given element or compound of a given stoichiometry 
were shifted so that to match the experimental cohesive energy of the equilibrium ground 
state structure. This procedure is followed merely to facilitate the comparison of structural 
energies calculated by different methods and does not introduce any new approximation. 
The equilibrium energy of the DOig-TiaAl compound has also been calculated here to check 
the methodology and provide useful reference information. The energy of this phase was 
minimized with respect to its volume while keeping the c/a ratio fixed at the value found in 
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The EV curves were calculated using the full potential LAPW metho d 46l47l4S within the 
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Hohenberg-Kohn-Sham formulation of the density functional theory . 49 i 50 i 51 The calculations 
were carried out in a spin-restricted mode and the exchange-correlation effects were treated 
at the level of the GGA. The EV calculations for Ti were carried out using an LAPW code 
arable at the Nava. Research Laboratory and were reported hr Ref.Q The calculations for 
Al, TiAl and Ti 3 Al were performed in this work and employed the WIEN2K package. 52 The 
Perdew-Wan g 53 i 54 i 55 (PW91) exchange-correlation functional was used for Ti calculations,— 
while its simplified and more efficient version referred to in literature as PBE^ was used for 
Al and TiAl. The muffin-tin radii of Al and Ti were chosen to be 2 a.u. For each crystal 
structure, systematic k point and basis set convergence tests were carried out at a fixed 
volume near equilibrium volume. The same set of parameters was subsequently employed 
for different volumes as well as for the Bain path calculations. The accuracy of the LAPW 
total energies calculated in the present work was estimated to be better than 0.5 mRy/atom. 

We now describe some relevant details of the potential fitting procedure. In the embedded 
atom formalism , 17 ! 18 the total energy of a system is expressed as 

Here, <3>jj is the pair- interaction energy between atom i and j at positions fi and r}, and Fi 
is the embedding energy of atom i. The pi in Eq. is the host electron density at site i 
induced by all other atoms in the system. The latter is given by 

Pi = J2Pj( r ij)- (2) 
For a binary system A-B, the total energy given by Eq. (0) is invariant with respect to the 
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following set of transformations : 31 ' 

p A (r) -> s A p A (r), (3) 

Pb(t) -> s B p B (r), (4) 

- F A {p A (r)ls A ), (5) 

Fs(p) - F B (p B (r)/s B ), (6) 

F A (p) ^ F A (p) + g A p, (7) 

F B (p)^F B (p)+g B p, (8) 

$A4(r) -> $A^(r) - 2g A p A (r), (9) 

-> - 2<7 B p B (r), (10) 

-> $ AB (r) - 2g A p A (r) - 2g B p B (r), (11) 



where A and S refer to the type of element (Al or Ti) and s A , s B , g A , and g B are arbitrary 
constants. For any particular compound, the potential can be cast into the so-called effective 
pair format 57 by choosing g A = —F' A (p~ A °) and g B = —F' b (pb ), where p A ° and p~ B ° are the 
equilibrium electron densities on atoms in the compound. The effective pair format provides 
a convenient way of comparing different potentials for the same compound. 

An EAM potential for a binary system can be constructed by two different procedures. 
One is to optimize all potential functions simultaneously during a single fit, as it was done for 
example for the Ni-Al system.— This scheme offers the advantage of having many parameters 
available for fitting to properties of a selected alloy or compound. However, it suffers from 
the drawback that the quality of potentials for pure elements (Al and Ti in our case) is often 
poor. An alternative and more common approach is to separately develop accurate EAM 
potentials for pure elements and use them to fit the cross potential for alloys (compounds). 
We resort to the latter procedure for the Ti-Al system. The parametrization of potential 
functions employed in the fitting procedure is discussed below. 

B. EAM potential for Al 

For the EAM potential for Al, we chose the electron density function in the form 

p(r) = i/j {^p 1 ^ {Mr ~ r fe-^-^ 

x (l + J B oe -^-^)+Co}, (12) 
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Here, A Q , B , Co, r , r c , h, y and 7 are the fitting parameters and ip(x) is a cutoff function 
defined as 



ip(x) = for x > 



x 



1 + x 4 



for x < 0, 



(13) 



where r c is the cutoff distance. The electron density in equilibrium fee Al is normalized to 
unity, i.e. p = ^\ NjPj = 1 , where j runs over coordination shells and Nj is the number of 
atoms on the j-th coordination shell. This constraint fixes one parameter in the above set 
of parameters. The pair interaction function is parametrized in the form 



V 



(62-61) \z 



.61 



z b 2 



+ 5 



4> 



r — r, 



h 



(14) 



where, z = r/r', and 61, 62, 5, Vq, and r' are fitting parameters. Thus, altogether we have 
12 fitting parameters for functions p(r) and $(r). 

The embedding energy is obtained by equating the energy of fee Al [Eq. ((TJ)] to the 
universal equation of state (EOS). By studying a broad range of materials, Rose et alM 
proposed the "universal" EOS in the form 



E(r) = -E [l + a(r/r e - l)]e 



-a(r/r e — 1) 



(15) 



where a = ^QVLqB / Eq, and r, Eq, B, r e are the nearest-neighbor distance, cohesive 
energy, equilibrium atomic volume, the bulk modulus, and the equilibrium-nearest neighbor 
distance, respectively. It is generally found that the EAM potentials which exactly follow 
Rose's EOS [Eq. ((15])] underestimate energies at high pressures.— We have therefore modified 
Eq. (JT5|) to allow for a more accurate fit to energies at high pressures. The modified equation 
has the form 



E(r) 



-E 



1 + ax + (5a 6 x 6 



(16) 



(x + 1) 2 

with x — (r/r e — 1). The parameter (3 in this equation is related to the pressure derivative 
of the bulk modulus at equilibrium as B' = |a + Qotf3 + 2. This modification does not alter 
the exact fitting of the potential to Eo,r e and B but provides a way to achieve an accurate 
fit to the experimental pressure-volume relation by adjusting the value of (3. 



7 



C. EAM potential for Ti 



Titanium has the hep structure at T = (a-Ti phase). For hep metals, an EAM potential 
can only be fitted to elastic constants if the relation (3Ci 2 — Cn)/2 > (Ci 3 — C 44 ) is 
satisfied.— Fortunately, this relation holds for a-Ti. 

In the present EAM potential for Ti, the electron density function is described by 



p{r) 



y[ e -ai{r-ro) 2 _|_ e ~a2(r-r' ) 



■0 



r — r r 



h 



(17) 



where the cutoff function ip(x) is given by Eq. ()13|) with the fitting parameters 
A,ai,a2,r ,r' ,r c , and h. One of the parameters is fixed by the normalization condition 
p — 1 at equilibrium. 

The pair interaction function is represented by 



r — r. 



h 



Voe -Ptr 1+V , 



-2l3 2 (r-r[) _2 e ^{r-r' 1 ) 



+ 8>, 



(18) 



where Vb, Vq, Pi, fa, r i, r [, and 5 are fitting parameters. The embedding energy function is 
expressed ClS db polynomial: 



F(p) = F + -F 2 (p-l) 2 + g (p-l) 3 



i+3 



(19) 



i=i 



Here F and F 2 are the embedding energy and its second derivative at equilibrium, respec- 
tively. These can be expressed in terms of the experimental values of E , B, and Q Q : 



F Q = E -±J2 N & 



and 

3 3 

where j runs over coordination shells, Nj is the number of atoms on the j-th coordination 
shell of radius Rj, while $j and <f>" are the pair interaction energy and its second derivative 
evaluated at Rj. The coefficients go and E>i {i < 3) in Eq. (pP?]) are fitting parameters. The 
parameter q was adjusted to ensure that the embedding energy vanishes when the electron 
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density goes to zero, that is, F(0) = 0. This requirement leads to the following expression 
for q : 

q = F + ^ + B 1 -B 2 + B 3 . (20) 

During the optimization of fitting parameters, the energy of the hep structure was required 
to approximately follow the Rose's EOS [Eq. ()15j) ] in the neighborhood of equilibrium. This 
was achieved by adding to the objective function the mean-squared deviation of the energy 
from Eq. (fTo^) at several points near the equilibrium. 

D. The cross potential Ti-Al and the fitting procedure 

Once the EAM potentials for Al and Ti are obtained, the cross potential representing the 
interactions between Ti and Al atoms was constructed by employing the parametrization 
given by Eq. (JTJJ). The transformation coefficients SAi,gn, and g^i [see Eqs. (3)-(ll)] were 
used as additional adjustable parameters. 

It should be mentioned that the specific analytical forms of the potential functions 
adopted in this work were found by trying a number of different forms and selecting those 
which provided a better accuracy of fitting with less parameters. The optimized values of 
the fitting parameters are listed in Table HI The potential functions are plotted in Fig. Q 
in the effective pair forma t? 1 ^ with respect to 7-TiAl . These functions are available in 
the tabulated form on Internet^ or from the authors upon request. The cutoff radii of 
atomic interactions in Al, Ti, and TiAl are 6.72, 5.19, and 5.77 A, respectively. The fitting 
procedure involves the total of 39 independent fitting parameters. 

III. TESTS AND APPLICATIONS OF THE POTENTIALS 
A. Aluminum 

The accuracy of the fitted EAM potential for Al can be adjudged from Table ITT] wherein 
the basic lattice properties, elastic constants, vacancy formation and migration energies, 
surface energies, and the stacking fault energy predicted by the potential are compared with 
their experimental values. We have also included the results obtained with our previous EAM 
potential for Al, 28 which we hereafter refer to as MFMP. The results obtained with the new 

9 



EAM potential are in good agreement with their experimental counterparts. The calculated 
vacancy formation energy El and migration energy E™, which are important for studying 
point defect diffusion, are well reproduced. The calculated vacancy formation volume fl( 
compares well with the one obtained with the MFMP potential as well as with experimental 
data. 60 The predicted intrinsic stacking fault ('Jsf) energy is on the lower end of the range 
of experimental values while the ab initi o 101 ' ' 102 calculatations yeilds a higher value of 166 
mJ/m 2 . The calculated symmetrical twin boundary (77-) energy is in good agreement with 
its experimental counterpart. The unstable stacking fault energy is underestimated with 
respect to the ab initi o 101 ' ' 102 value of 220 mJ/m 2 . The surface energies are underestimated 
in comparison with experiment, which is a general characteristic of EAM potentials. 

The structural energy differences for Al, given in Table 11111 reproduce ab-initio energies 
reasonably accurately. The EAM potential predicts the c/a ratio of the hep structure to be 
1.63 in good agreement with the optimized value of 1.645 obtained by the LAPW calcula- 
tions. In agreement with earlier findings,— bec Al is mechanically unstable (Cn < C12) and 
transforms to the fee structure upon c/a relaxation. The sc structure is also mechanically 
unstable with C44 < 0. A comparison of the EOS's of various crystalline structures of Al 
calculated with the EAM potential and by the LAPW method is presented in Fig. |21 The 
EAM curves are seen to agree with the LAPW results fairly well. The agreement is particu- 
larly good for the bec and fee structures over a large range of volumes, but tends to worsen 
for more open structures. 

As was mentioned in Section II B, the Al potential was fit to the experimental P(V) 
relation by adjusting the parameter f3 in the generalized EOS, Eq. ()16j) . The optimized value 
of P = 0.00489 provides an excellent agreement with experimental data up to pressures of 
about 700 GPa, as illustrated in Fig. 01 In contrast, the standard Rose's EOS (/3 = 0) 
underestimates the pressures under strong compressions. Note that both equations share 
the same values of Eq, B and Oq. The excellent fit to high pressure data makes the potential 
useful for simulating shock waves, sputter deposition, and other processes involving a close 
approach of atoms. 

Thermal expansion of Al was studied within the temperature range of 5-1000 K. The 
calculated thermal expansion factors at selected temperatures are given in Table IIV1 They 
were obtained using a 864-atom supercell by two different methods. In the first method, 
the free energy of the crystal was minimized as a function of volume in the quasiharmonic 
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approximation.— This method includes quantum-mechanical effects such as zero-point vi- 
brations and should yield more accurate values of the thermal expansion factor at low tem- 
peratures. However, it may not be very accurate at high temperature where the anharmonic 
effects become significant. The second type of calculation was carried out by the Metropolis 
Monte Carlo method . 61 ' 62 This method is based on classical mechanics and fully incorporates 
anharmonic effects. It is therefore more adequate for thermal expansion calculations at high 
temperatures. As can be seen from Table ITVl the MC results are quite close to experimental 
data at high temperatures. 

Overall, the EAM potential developed here provides a good description of a wide range 
of Al properties. Despite the existence of other high-quality EAM potentials for Al in the 
literature , 23 i 28 i 34 we chose not to re-use one of them but rather generate a new potential so 
that to have all potential functions for the Ti-Al system created by the same methodology. 
We also used this work as an opportunity to address some weak points of previous potentials. 
For example, even though the MFMP potential demonstrates a better agreement with 
experiment for some of the properties listed in Table |H1 the present potential describes the 
thermal expansion and high-pressure behavior of Al more accurately. The present potential is 
also based on a larger set of ab-initio data and should be better transferable to configurations 
away from equilibrium. It should also be mentioned that the use of smooth analytical 
functions in this work makes the potential more robust in comparison with the cubic-spline 
parameterization applied in Ref. I28L 



B. Titanium 



Equilibrium lattice properties, vacancy characteristics, as well as stacking fault and sur- 
face energies in Ti computed using the present EAM potential are compared with experi- 
mental data in Table IVl We have also included the results obtained with the EAM potential 
developed by Fernandez et al. , 37 which will be referred to as FMP. The latter is an improved 
version of the potential proposed in Ref. |35j. Although our potential was not fitted to the 
c/a ratio exactly, the predicted value of 1.585 is in a good agreement with the experimental 
value of 1.588. The elastic constants are also reproduced reasonably well. 

The vacancy formation energy El was fitted to the target value of 1.85 eV. The exper- 
imental value of El reported by Shestopal 63 is 1.55 eV, more recent positron annihilations 
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measurements^ give E{ = 1.27 eV, while the ab-initio linearized muffin-tin orbital method 
(LMTO) method 6 ^ yields a much higher value of 2.14 eV. We therefore opted to fit to an 
intermediate value of 1.85 eV, which after the relaxation decreased to 1.83 eV. The vacancy 
migration energy E™ was calculated using the nudged elastic band method.— In the hep 
lattice with a nonideal c/a ratio, the basal and non-basal vacancy jumps are not identical. 
The calculated values of E™ for the basal and non-basal jumps are 0.80 eV and 0.83 eV, 
respectively. The FMP potential gives smaller values of 0.51 eV and 0.48 eV, respectively. 
While the experimental value of the vacancy formation energy is rather uncertain, the ex- 
perimental activation energy Q of self-diffusion in a-Ti , which is the sum of El and E™ 
, has been measured fairly accurately 67 For self-diffusion perpendicular to the c axis, the 
experimental value is Q = 3.14 eV. The present EAM potential predicts Q = 2.62 eV, while 
the FMP potential gives a lower value of 2.02 eV. 

There are three stacking faults on the basal plane in a-Ti , which are deviations from the 
normal stacking sequence ABABAB of closed packed planes in the hep structure . 68 ' 69 The 
intrinsic fault I\ is formed by a removal of one hexagonal layer followed by a | < 1010 > 
translation of all atoms above this fault. The resultant stacking sequence is 

h : ABAB\CBCBC, 

where the vertical bar indicates the position of the fault. The intrinsic stacking fault I2 is 
created by a | < lOlO > slip: 

I 2 : ABAB\CACAC. 

The extrinsic stacking fault E result from the insertion of an extra hexagonal plane into the 
normal stacking sequence: 

I E : ABAB\C\ABAB. 

The calculated relaxed stacking fault energies (Table EJ compare well with those obtained 
with the FMP potential. The experimental value of the I2 fault energy is about 300 mJ/m 2 
and is considered to be a rough estimate. Both EAM potentials underestimate this exper- 
imental value. All our efforts to obtain a higher 7j 2 value during the fitting of the present 
potential did not have much success. In fact, any attempt to raise ji 2 above 66 mJ/m 2 
resulted in a deterioration of other properties, which gave us an indication that higher 
stacking fault energies may be beyond the capabilities of the EAM. Note, however, that the 
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EAM-predicted stacking fault energies follow the expected relation^ Je ~ §7/ 2 ~ 3jh- The 
limited success in fitting to higher stacking fault energies is likely to be due to the directional 
component of bonding in Ti owing to d electrons. The covalent nature of bonding cannot be 
described by the central-force-based EAM model. More rigorous parameter-based methods 
such as the modified EAM ffiZS bond order potential o? 1 ^ 2 ^? or the tight-binding method/ 4 
which include angular-dependent interactions, may give higher stacking fault energies. 

The predicted value of the (0001) surface energy, 1725 mJ/m 2 , slightly underestimates 
the experimental value. This is again consistent with the general trend of the EAM to 
underestimate surface energies. The FMP potential yields an even smaller value of 1439 
mJ/m 2 , whereas the ab-initio surface energies, 2100 mJ/m 2 (Ref. [3) and 1920 mJ/m 2 
(Ref. 76), overestimate the experimental value. 

The LAP W and the EAM energies of various crystal structures of Ti relative to the hep 
structure are reported in Table IVI1 We note that the LAPW calculations predict the omega 
structure to be the ground state, with the hep energy being 0.06 eV/atom higher. We have, 
therefore, excluded the omega structure from the fitting procedure. Overall, both EAM 
potentials yield similar energy differences between the structures, with the present potential 
performing somewhat better. Both potentials predict the hep structure to be more stable 
than the omega structure in agreement with experiment. In Fig. the EOS's of the hep, 
fee, bec, and the sc structures of Ti calculated with the present EAM potential are compared 
with the LAPW results. The agreement between the two calculation methods is good for 
the close-packed structures but becomes poorer for the low coordinated sc structure. In the 
latter case, however, the present EAM potential is closer to the LAPW data than the FMP 
potential. 

The linear thermal expansion of Ti was calculated within the quasiharmonic approxima- 
tion and by the MC method using a supercell with 800 atoms. The c/a ratio was kept fixed 
at its equilibrium T = value during the calculations. The obtained values of the ther- 
mal expansion factor for selected temperatures are reported in Table IVHI The agreement 
with experimental data^ for polycrystalline Ti is reasonable. The FMP potential gives a 
poorer agreement with experiment. For example, at 293 K the FMP potential gives the 
quasiharmonic linear thermal expansion of 1.35% while the experimental value is 0.15%. 
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C. Intermetallic compound 7-TiAl 



The physical properties of 7-TiAl obtained with the present EAM potential are summa- 
rized in Table IVTTT1 The lattice constant and the cohesive energy are reproduced accurately. 
The cj a ratio is correctly predicted to be larger than unity and is in good agreement with the 
experimental value. The elastic constants are in a reasonable agreement with experiment, 
the root-mean-squared deviation for elastic constants being about 22%. We note that the 
negative signs of the two Cauchy pressures, (C12 — Cqq) and (C13 — C44), are not reproduced 
by the present potential, nor are they reproduced by previous EAM-type potentials. The 
negative Cauchy pressures in TiAl are caused by the directional component of bonding and 
cannot be described by the EAM. Table lyiTTl also includes the numbers calculated with the 
Farkas potentia l?^ as well as with the P2 potential constructed by Simmons et al. 21 The 
respective root-mean-squared deviations of the elastic constants from experimental data are 
45% and 28%. 

The planar defect energies in 7-TiAl are summarized in Table IIXI They were calculated 

1 /3 

using supercells with an effective cubic lattice with the lattice parameter a = (a^Co) an d 
without volume relaxation. Since the cohesive energy of the equilibrium tetragonal lattice 
was used as a reference in the calculations, the resultant fault energies can be slightly over- 
estimated. For the superlattice intrinsic stacking fault (SISF), the EAM value of 173 mJ/m 2 
is slightly higher than the experimental value of 140 mJ/m 2 . There have been a number 
of ab-initio calculations of the SISF energy with results scattered over the range of 90-172 
mJ/m 2 . For the complex stacking fault (CSF) energy we obtain the value of 299 mJ/m 2 
well bracketed between ab-initio results. Experimental data for the CSF energy are not 
available. The calculated value of the antiphase boundary (APB) energy, 266 mJ/m 2 , is 
also in a good agreement with the experimental value of 250 mJ/m 2 . The ab-initio APB 
energies are scattered over the wide range 510-670 mJ/m 2 . The hierarchy of planar fault 
energies in 7-TiAl was investigated by Wiezorek and Humphreys.-^ According to their pre- 
liminary computational results, this hierarchy in Ti-54at%Al is ■Jcsf > Iapb > Isisf- The 
present EAM potential predicts the same ordering. On the other hand, ab-initio calcula- 
tions with the linearized Korringa-Kohn-Rostoker (LKKR) method 80 and the full potential 
LAPW (FLAPW) method 81 give the japb > Icsf > Isisf ordering. More recent calcula- 
tions by Ehmann and Fahnle by the LAPW method including local atomic relaxations are 
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consistent with the latter ordering of the stacking fault energies. The discrepancy between 
the experimental data and EAM calculations, on one hand, and ab-initio calculations, on 
the other hand, originates primarily from the high APB energy delivered consistently by 
ab-initio methods. The low APB energy observed experimentally may reflect the local dis- 
order near the APB taking place due to the off-stoichiometry and/or temperature effects. 
The similarly low APB energy predicted by the present EAM potential can lead to a good 
agreement between atomistic simulations and experiment. 

Simmons et ah — succeeded in generating a set of EAM potentials for 7-TiAl fit to high 
APB energies comparable to ab-initio values, but their potentials give c/a < 1 in contradic- 
tion to experimental data. When generating our potential we could also achieve higher ^apb 
values at the expense of c/a < 1, but could never increase Japb above 266 mJ/m 2 while 
keeping c/a > 1 and maintaining a good quality of fit to other properties. We believe that 
the underestimation of the APB energy is another intrinsic limitation of the central-force 
EAM as applied to 7-TiAl . Farkas^ constructed a potential that gives japb > 1c sf while 
c/a > 1. However, some of the elastic constants predicted by that potential are in a poor 
agreement with experimental data, especially C13 and C33 (cf. Table IVrTT|) . The potential 
also gives a discontinuous temperature dependence of the quasiharmonic thermal expansion, 
with unrealistically large values at high temperatures. 

The EV curves computed with the present EAM potential and by the LAPW method are 
presented in Fig. El The agreement between the two calculation methods is good for the Llo 
structure. The difference between the two curves for the B2 and Bl structures is presumably 
due to the limited accuracy of the EAM method to describe the open structures. The 
formation energies for different structures of TiAl (relative to fee Al and hep Ti) calculated 
by the LAPW method and with the present EAM potential are presented in Table [Xj For 
comparison, experimental and ab-initio results available in the literature have also been 
included in the Table. The formation energies obtained with the present EAM potential 
are in good agreement with the corresponding experimental and ab-initio energies. The 110 
structure is correctly produced to be the ground state. We note the B2 and 532 structures 
are unstable with respect to the c/a optimization. In particular, the B2 structure transforms 
to the equilibrium 110 phase upon c/a relaxation. 

For Ti 3 Al , the EAM potential correctly predicts the DO19 structure to be the equilib- 
rium ground state of Ti 3 Al (Table EJ). The lattice constants, cohesive energy, and the 
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elastic constants of DO19-T13AI are given in Table IXII We emphasize that none of these 
properties were included in the potential fit. The observed agreement with experimental 
data demonstrates a good transferability of our potential. 

The EAM potential was also applied to investigate the Al 3 Ti compound. Experimentally, 
the equilibrium structure of Al 3 Ti is DC>22- The present EAM potential predicts the Ll 2 
structure to be 0.01 eV lower in energy than the DO22 structure, suggesting that the potential 
may not suitable for simulating the Ah>Ti compound. 

Thermal expansion factors of 7-TiAl calculated within the quasiharmonic approximation 

able IXHI The calculated values are in 



82 



and by the Monte Carlo method are presented in 
agreement with those estimated from Fig. 4 in Ref. 

The energy along the Bain path between the tetragonal 7-TiAl structure and the B2 
structure was calculated by the EAM and LAPW methods. Starting from the equilibrium 
tetragonal 7-TiAl structure, the c/a ratio was varied by keeping the volume constant. The 
energy change during the transformation is plotted in Fig.|H]as a function of the deformation 
parameter X defined by c/a = X(co/ao)- The EAM energies are observed to closely follow 
the LAPW energies along the path. This agreement confirms a good transferability of the 
present EAM potential. 

Point defect properties play an important role in the atomic disorder and diffusion in 
7-TiAl . The TiAl lattice supports two types of vacancy (Vti and Vai) and two types 
of antisite defects (Ti atom on the Al sublattice, Tiai, and Al atom on the Ti sublattice, 
Afpi). 12 The so-called "raw" formation energies and entropies^ 2 - of the defect formation have 
been calculated with the present EAM potential using the molecular statics method for the 
energies and the quasiharmonic approximation for the entropies. 

When analyzing point defects in ordered compounds it is more convenient to deal with hy- 
pothetical composition-conserving defect complexes rather than individual defects . 12 i 83 i 84 i 85 
It should be emphasized that the defects are grouped into complexes conceptually and not 
physically. The complexes are assumed to be totally dissociated and interactions between 
their constituents are neglected. The advantage of dealing with composition-conserving 
complexes is that all reference constants involved in their energies and entropies cancel out. 
This allows us to directly compare results obtained by different calculation methods. The 
complex energies and entropies can be expressed in terms of the "raw" energies and en- 
tropies Sd, d = Vxi, Vai, Tiai, Al Ti . 12 The expressions for some of the complex energies are 
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given in Table IXIII1 Similar expressions hold for the complex entropies, except that the 
cohesive energy E Q should be replaced by the perfect lattice entropy per atom. Table IXlTTl 
summarizes the results of the EAM calculations for several defect complexes and compares 
them with the ab-initio energies reported by Woodward et a/.— The agreement between the 
two calculation methods is reasonable. We emphasize again that point defect properties of 
7-TiAl were not included in the potential fit. 

Using the complex energies and entropies, the equilibrium defect concentrations have 
been calculated as functions of the bulk composition around the stoichiometry within the 
lattice gas model of non-interacting defects . 12 i 83 i 84 i 85 Fig. shows the calculation results 
for T = 1000 K. We see that all compositions are strongly dominated by antisite defects. 
This observation is well consistent with the experimentally established fact that 7-TiAl 
is an antisite disorder compound . 87 ' 88 ' 89 The vacancy concentrations are several orders of 
magnitude smaller than antisite concentrations. In the stoichiometric composition, most 
of the vacancies reside on the Ti sublattice. All these features have been observed at all 
temperatures in the range 800-1200 K. 

IV. SUMMARY 

EAM potentials have been developed for Al, a-Ti and 7-TiAl by fitting to a database of 
experimental data and ab-initio calculations. The potentials have been tested against other 
experimental and ab-initio data not included in the fitting database. The ab-initio structural 
energies for Ti were calculated previously— while those for Al and Ti-Al compounds have 
been generated in this work. All these calculations employed the full-potential LAPW 
method within the GGA approximation. Besides serving for the development of the EAM 
potentials, the obtained ab-initio energies are also useful as reference data for the Ti-Al 
system. 

The Al potential is fit to the target properties very accurately and has demonstrated 
a good performance in the tests. It has certain advantages over the previously developed 
potential,— particularly with respect to the lattice thermal expansion and the pressure- 
volume relation under large compressions. The fit of the Ti potential is less successful, 
presumably because of the directional component of interatomic bonding that is not captured 
by the central-force EAM model. In particular, the potential underestimates the stacking 
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fault energies on the basal plane. Further improvements of the potential do not appear to be 
possible within the EAM. This potential can be viewed as a supporting potential for the Ti- 
Al system, but we also believe that it can be useful in atomistic simulations in pure Ti where 
subtle details of atomic interactions may not be critical. Since the potential is fit reasonably 
well to the elastic constants, thermal expansion factors and the vacancy formation energy, 
it can be employed for modeling diffusion and creep in large systems that are not accessible 
by more accurate, yet slower, ab-initio methods. 

For the 7-TiAl compound, the potential developed here reproduces reasonably well the 
basic lattice properties, planar fault energies, as well as point defect characteristics. The 
fit to the elastic constants is better than with previous potentials. However, the negative 
Cauchy pressures in 7-TiAl have not been reproduced by the present nor previous EAM 
potentials. The planar fault energies calculated with the potential are in a good agreement 
with experiment, but the APB energy is lower than all ab-initio values. The fit to ab-initio 
energies of alternative structures of TiAl enhances the transferability of the potential to 
configurations away from equilibrium. This fact is verified by the good agreement between 
the EAM and LAPW energies along the Bain transformation path. The potential also cor- 
rectly predicts the equilibrium DO19 structure of Ti 3 Al and gives a fairly good agreement 
with experiment for the cohesive energy, lattice parameters, and elastic constants of this 
compound. The point defect energies and entropies in 7-TiAl calculated with the potential 
are in agreement with the antisite disorder mechanism established for this compound ex- 
perimentally. We emphasize that neither the Bain path nor any information on TisAl or 
point-defect properties in 7-TiAl were included in the fitting database. This success of the 
proposed potential points to its ability to describe atomic interactions in the Ti-Al system 
on a reasonable quantitative level. The potential should be suitable for large-scale atomistic 
simulations of plastic deformation, fracture, diffusion, and other processes in 7-TiAl . At the 
same time we acknowledge that more rigorous models, particularly those including angular- 
dependent interactions, are needed for addressing the negative Cauchy pressures, high APB 
energy, and other properties of 7-TiAl that lie beyond the capabilities of the EAM. 
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FIG. 1: The embedding energy (a), electron density (in arbitrary units) (b), and the pair interac- 
tion function (c) for the Ti-Al system in the effective pair format. 
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FIG. 2: Energy- volume relations for different crystalline structures of Al calculated with the EAM 
potential (lines) and by the LAPW method (points). 
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FIG. 3: The pressure- volume relation for aluminum at T = calculated with the present EAM 
potential (solid line), predicted by the universal equation of stated (dotted line), and measured 



experimentally (squares: Ref. l90l : circles: Ref. l9ll : triangles: Ref. 
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FIG. 4: Comparison of the LAPW (points) and EAM (solid lines) energy-volume relations for 
different structures of Ti. The LAPW energies were calculated in Ref. |4J. 
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FIG. 5: Energy- volume relations for the LIq, B2 and Bl structures of TiAl. The LAPW results 
are marked by points, the lines represent EAM calculations. 
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FIG. 6: Energy per atom as a function of the deformation parameter X (see text for details) 
along the volume conserving tetragonal deformation path (Bain path) in TiAl. The energy is 
given relative to the equilibrium LIq structure. The dotted line is predicted by the present EAM 
potential and the symbols represents the LAPW results. 
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FIG. 7: Calculated equilibrium concentrations of vacancies and antisites in TiAl as functions of 
the alloy composition at 1000 K. 
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TABLE I: Optimized values of the fitting parameters of the EAM potential for the Ti-Al system. 
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r' ( k\ 


9 ^Qfi^A^ v in 2 








0.00489 


«i (i/A) 


0.3969775 x lO" 1 










a 2 (1/A) 


5.344506 x 10 2 










Bi 


1.549707 










B 2 


-0.4471131 










B 3 


0.8594003 x lO" 1 
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TABLE II: Properties of Al calculated using the present EAM potential and the MFMP potentia 
in comparison with experimental data. 



28 



Property 



Experiment 



EAM 



MFMP 



Lattice properties: 

a (A) 

Eo (eV/atom) 

B (GPa) 

Cu(GPa) 

Ci 2 (GPa) 

C 44 (GPa) 

Vacancy: 

E^(eV) 

E-(eV) 

nl/n Q 

Planar defects: 
75F(mJ/m 2 ) 
7 us (mJ/m 2 ) 
7T(mJ/m 2 ) 
Surface: 

7 s (110)(mJ/m 2 ) 
7 s (110)(mJ/m 2 ) 
7 s (lll)(mJ/m 2 ) 



4.05 a 
3.36 fe 

79 c 

114 c 
61.9 C 
31.6 C 

0.68^ 
0.65 e 
0.62^ 

166 s ,120-144' 1 

769 

980*, lUW, 1160 fc 
980*, 1140- 7 , 1160 fc 
980*, 1140*', 1160 fc 



4.05 
3.36 
79 

116.8 

60.1 

31.7 

0.71 
0.65 
0.59 

115 
151 
63 

792 
607 
601 



4.05 
3.36 
79 

113.8 

61.6 

31.6 

0.68 
0.64 
0.51 

146 
168 
76 

1006 

943 

870 



b Ref. m 
c Ref. m 
d Ref. m 
e Ref. m 
^Ref. m 
s Ref. M 
h ReL ,99 and 
Average orientation, Ref. 98 
^Average orientation, Ref. ft 
k Average orientation, Ref. 1 75 



i lionl 
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TABLE III: Comparison of the energies (eV/atom) of selected structures of Al calculated by the 
LAPW method, with the present EAM potential, and with the MFMP potential. — The energies 
are given relative to the energy of the equilibrium fee structure. 

Structure LAPW EAM MFMP 

hep 0.04 0.03 0.03 

bee 0.09 0.09 0.011 

Ll 0.27 0.33 0.30 

sc 0.36 0.30 0.40 

diamond 0.75 0.88 0.89 
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TABLE IV: The linear thermal expansion factor (in %) of Al computed using the present EAM 
potential in comparison with experimental data at selected temperatures. QHA: quasiharmonic 
approximation; MC: Monte Carlo method. 



T(K) 


Experiment 


QHA 


EAM 

MC 


293 


0.418 


0.277 


0.489 


500 


0.932 


0.663 


0.872 


700 


1.502 


1.016 


1.332 


900 p- . 


2.182 


1.419 


1.916 


"Ref. 77 
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TABLE V: Properties of Ti predicted by the present EAM potential and the FMP— potential in 
comparison with experimental data. 







EAM 


FMP a 


an (A) 


2 951 a 


2 951 


2 951 


c/a 


1 5R8 a 

1 . (JOO 


1 585 


1 588 


F,^ ('pV/fltnml 

HlQ tc V / dLLUil 1 


4 850 fe 

4:.OtJU 


4 850 


4 850 




1 7fi c 


1 78 

1 1 o 


1 8Q 




87 c 


74 


74 




uo 


77 


UO 


^33 


1 QD C 


1 Q1 


xoo 


C^TOPa! 




51 


5(1 




1 


1 8^ 


1 51 


E m fhasal") (pV) 




n so 


n 51 


E 1 ™ (nonbasal) (eV) 




0.83 


0.48 


Q (eV) 


3.14^ 


2.62 


2.02 


7/j (mJ/m 2 ) 




31 


31 


7/ 2 (mJ/m 2 ) 


290 5 ,300 /l 


56 


57 


(mJ/m 2 ) 




82 


84 


7s (0001) (mJ/m 2 ) 


2100*, 1920?' 


1725 


1439 



°Ref. Wfl 
6 Ref. 22 
c Ref. gS 
d Ref. S3 
^Ref. E 
»Ref. 
h Ret. 



104 



105 



* Average orientation, Ref. 75 
^Average orientation, Ref. 176 
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TABLE VI: Energies (eV/atom) of selected structures of Ti obtained with the present EAM 

potential and by the LAPW/GGA-PW91 method^ All energies are given relative to the energy 
of the experimentally observed hep structure. 

Structure LAPW/GGA-PW91 EAM FMP 

fee 0.012 0.011 0.012 

bec 0.067 0.03 0.02 

sc 0.77 0.54 0.27 

omega -0.06 0.094 0.064 
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TABLE VII: The linear thermal expansion factor (in %) of Ti calculated with the present EAM 
potential in comparison with experimental data at selected temperatures. QHA: quasiharmonic 
approximation; MC: Monte Carlo method. 



T(K) 


Experiment" 


EAM 

QHA 


MC 


293 


0.15 


0.16 


0.25 


500 


0.35 


0.38 


0.44 


700 


0.55 


0.62 


0.63 


100Q— . 


0.89 


1.00 


0.72 



"KFJTTT 
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TABLE VIII: Equilibrium lattice constant, c/a ratio, cohesive energy, and elastic constants of 
7-TiAl calculated with the present EAM potential in comparison with other potentials 21,38 and 
experimental data. 



Experiment EAM Ref.38 Ref.Jl 



a (A) 




3.997 a 


3.998 


3.951 


4.033 


c/a 






1.02 a 


1.047 


1.018 


0.991 


Eq (eV/atom) 


4.51 6 


4.509 


4.396 


4.870 


Cu(Gp a ; 




186 c , 183 d 


195 


222 


202 


Ci 2 (Gp a ; 




72 c , 74.1 d 


107 


100 


95 


Ci 3 (Gp a ; 




74 c , 74A d 


113 


162 


124 


C 33 (GPa; 




176 c , 178 d 


213 


310 


237 


C 44 (GPa; 




101 c , 105 d 


92 


139 


83 


C 66 (G 






77 c , 78.4 d 


84 


76 


54 


"Ref. 
b Rei. 
c Ref. 
d Ref. 


100 

m 

108 
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TABLE IX: The energies of the superlattice intrinsic stacking fault (SISF), antiphase boundary 
(APB), complex stacking fault (CSF), and the (100) and (111) surfaces in 7-TiAl calculated with 
the present EAM potential. Available experimental and ab-initio data are included for comparison. 



All energies are expressed in mJ/m 2 . 

Experiment ab-initio EAM 

SISF(lll) 140° 90 6 , 110 a , 123 c , I72 d 173 

CSF(lll) 280 a , 294 c , 363 d 299 

APB(lll) 250 a 510 6 , 667 d , 670 a , 672 c 266 

surface (100) 1177 

surface (110) 1445 



^RefsfHy 
c Ref. tM 
d ReL 110 
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TABLE X: Formation energies, AH (eV/atom), of 
culated with the present EAM potential and by the 
data. 



different compounds of the Ti-Al system cal- 
LAPW method in comparison with literature 





TiAl 








T13AI 


Structure 


AH 


Method 


Structure 


AH 


Method 


Ll 


-0.404 


EAM a 


D0 19 


-0.289 


EAM a 


Ll 


-0.43 


LAPW/GGA a 


DO19 


-0.318 


LAPW/GGA" 


Ll 


-(0.37-0.39) Experiment 6 


D0 19 


-0.25,-0.26 


Experiment^ 


Ll 


-0.38 


Experiment 


DOi 9 


-0.28 


FLASTO/LDA d 


Ll 


-0.41 


FLASTO/LDA d 


DOi 9 


-0.29, -0.28 


LMTO/LDAf, FLAPW/LDA 3 


Ll 


-0.44 


FLMTO/LDA e 


DO19 


-0.28 


FLAPW/LDA' 1 


B 2 


-0.27 


EAM a 


L12 


-0.288 


EAM a 


B 2 


-0.29 


LAPW/GGA a 


L12 


-0.30 


LAPW/GGA a 


B 2 


-0.26 


FLASTO/LDA d 


L12 


-0.28 


FLASTO/LDA d 


Bi 


0.13 


EAM a 


L12 


-0.29, -0.27 


LMTO/LDA 5 , FLAPW/LDA 9 


Bi 


0.24 


LAPW/GGA a 


DOn 


-0.03 


EAM a 


Lla 


-0.30 


EAM a 


DO22 


-0.28 


EAM a 


B32 


-0.32 


EAM a 


DO22 


-0.27, -0.25 


LMTO/LDA 9 , FLAPW/LDA 3 


"40" 


-0.37 


EAM a 


D0 3 


-0.23 


EAM a 




Al 3 Ti 










Structure 


AH 


Method 








D0 22 


-0.29 


EAM a 








D0 22 


-0.41 


FLASTO/LDA d 








D0 22 


-0.42 


LMTO/LDA 8 








L12 


-0.30 


EAM a 








D0 3 


-0.20 


EAM a 








D0 19 


-0.29 


EAM a 









Treiprr 

6 Ref. 
c Ref. 
d Ref. 
e Ref. 



s 



wo 



113 
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^Refs 
sRef. n 
h Ref. tt 
'Ref. 15 
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TABLE XI: Equilibrium properties of T^Al predicted by the present EAM potential. Ah initio 
and experimental data are included for comparison. 



P T* C\ T\ V i~\ T 


iJAUCl IJ-LlCLiLi 


RAM 


nh—otiitin 

UiU Oil ILIU 


an ( k\ 

ao 


c 77a 


78.4 


^ fil 4 e 


c/a 




U.oZl 


U.ool 


Eo (eV / atom) 


4.78 6 


4.766 




Cu(GPa) 


176.2 C , 175 d 


180.5 


221/ 


Ci 2 (GPa) 


87.8 C , 88.7 d 


74.4 


71/ 


Ci 3 (GPa) 


61. 2 C , 62.3 d 


70.3 


85^ 


C 33 (GPa) 


218.7 C , 220 d 


222.9 


238^ 


C44(G£aj 


62.4 C , 62. 2 d 


46.6 


69^ 



6 Ref. 
c Ref. 



d Ref. rm 
e Ref. 4S 



/Ref. 



114 



107 



13 
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TABLE XII: The linear thermal expansion factor (in %) of 7-TiAl calculated with the present 
EAM potential. QHA: quasiharmonic approximation; MC: Monte Carlo method. 



T(K) 


QHA 


EAM 


MC 


400 


0.36 




0.53 


600 


0.66 




0.85 


800 


1.03 




1.20 


1000 


1.54 




1.58 
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TABLE XIII: Energies (in eV) and entropies (normalized to k^) of point defect complexes in TiAl 
computed using the present EAM potential. Ab-initio results are included for comparison. 



(~!oTTi"n1py 


Eoiiritinn 


nh-initin a 


- L -' iIci 6J 


EAM 


Divacancy 


ev Tl +ev Al +2E 


3.582 


3.168 


2.804 


Exchange 


eAhr t +eTi Al 


1.204 


0.765 


1.420 


Triple-Ti 


2ev Ti +2E Q +€Ti Al 


3.775 


3.525 


3.952 


Triple- Al 


2 e v Al +2E -\-eAi Ti 


4.593 


3.576 


3.075 


Interbranch Ti 


2 £v A r£Ti Al +2E 


3.389 


2.810 


1.656 


Interbranch Al 


eAi Ti -^ e v T i-^Eo 


-2.571 


-2.759 


-2.532 
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